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ABSTRACT 

In this model a collimated ultra-relativistic ejecta collides with an amor- 
phous dense cloud surrounding the central engine, producing gamma-rays via 
synchrotron process. The ejecta is taken as a standard candle, while assum- 
ing a gaussian distribution in thickness and density of the surrounding cloud. 
Due to the cloud high density, the synchrotron emission would be an instanta- 
neous phenomenon (fast cooling synchrotron radiation), so a GRB duration cor- 
responds to the time that the ejecta takes to pass through the cloud. Fitting the 
model with the observed bimodal distribution of GRBs' durations, the ejecta's 
initial Lorentz factor, and its initial opening angle are obtained as Fq < 10^, and 
Co ~ 10~^, and the mean density and mean thickness of the surrounding cloud 
as TT ~ 3 X lO^^cm"^ and L ~ 2 x lO^^cm. The clouds maybe interpreted as the 
extremely amorphous envelops of Thorne-Zytkow objects. In this model the two 
classes of long and short duration GRBs are explained in a unique frame. 

Subject headings: gamma rays: bursts — shock waves 



1. INTRODUCTION 

Undoubtedly, gamma-ray bursts (GRBs) have remained to be one of the most exit- 
ing, intriguing, and enigmatic astrophysical phenomena since their mysterious discovery in 
the past several decades (for a recent expository review of GRBs the interested reader is 
referred to the excellent work by J. I. Katz (2002)). Although no two GRBs resemble 
each other and each one has its own peculiarities which makes the problem of modeling 
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GRBs very difficult, the whole of GRBs reveals several interesting features. Since the pub- 
lication of the BATSE data (Fishman et al. 1993) which included the observation of over 
200 GRBs and revealed an almost uniform distribution of the location of GRBs in the 
sky, combined with the deficiency of faint GRBs, the association of GRBs with the galac- 
tic plane has been ruled out. The successive publications confirmed the figure more and 
more (Meegan et al. 1996; Paciesas et al. 1999). However, since the observation of after- 
glows in X-ray (Costa et al. 1997), optical (van Paradijs et al. 1997) and radio spectrum 
(Frail et al. 1997) and the advert of Robotic Optical Transient Search Experiment (ROT- 
SIE) telescope (Akerlof et al. 2003; Gisler et al. 1999, e.g.) which has revealed the red shifts 
for several GRBs, their cosmological origin is widely accepted. 

The BATSE data (Fishman et al. 1993: Meegan et al. 1996; Paciesas et al. 1999) has also 
revealed another equally important overall feature of the GRBs. The distribution of time 
duration in observed GRBs shows a double heap distribution, which the smaller one peaks 
around 0.2 sec and the larger one peaks around 20 sec (Fig.l). This two peaked distribution 
which apparently separated GRBs into the so called short duration and long duration ones 
was referred to as bimodality (Kouveliotou et al. 1993; Norris et al. 1993) and led some in- 
vestigators to beheve that there are two distinct populations of GRBs. 
It has been widely believed that whatever the central engine is, the radiation reaching us 
originates from the space surrounding the central engine. It is also believed that during the 
collapse the energy release streams out in relativistic confined ejectas (not isotropically) and 
thus the total energy release during each event is far less than the unbelievable amount that 
one might obtain by assuming isotropic radiation (Kulkarni et al. 1999). 
The aim of this paper is to present a rather simple model, based on these general ideas, and 
show that there is no need for assuming two distinct populations of GRBs, and to show that 
once the geometrical considerations and cosmological effects are fully accounted for a genetic 
standard candle ejecta, crossing an amorphous dense cloud, the so called bimodality can be 
deduced. 

In § 2 the model and its general formulation is introduced. In § 3 the computational results 
and the fit of model parameters with BATSE data will be presented as well. § 4 is devoted 
to a discussion on the results. Some needed calculations and discussions are presented in 
appendices. 

2. THE MODEL FORMULATION 

GRBs are modeled as a central engine with an instantaneous ultra-relativistic jet of 
material surrounded by an amorphous dense cloud. The central engine and its jet are taken 
as a standard candle with a total release energy of E, and an initial Lorentz factor Fq and 
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initial opening angle Co of the jet for all GRBs, whereas the clouds are considered to have a 
distribution both in thickness as well as in density. For the sake of illustration and brevity 
we take both of these distributions to be Gaussian. 

We want to calculate the distribution of logarithm of time duration of observed GRBs accord- 
ing to the above model. We should however explore the cjccta evolution since the observed 
gamma ray emission originates in the shock front of ejecta+shocked medium. 



2.1. The Ejecta Evoution 

The equations describing the ejecta evolution are presented here, based on the notations 
of Paczynski & Rhoads (1993). We consider the cloud to be at a distance Tq from the source, 
which may be negligible in comparison to the cloud thickness L. The ratio of swept-up mass 
to ejecta mass Mq has the form below: 

f = ^ f P^mir'y'dr', (1) 



where, r is the ejecta distance from the source. The cloud density p is taken to be independent 
of r. Furthermore, ^^^(r) = 27r[l — cosC(r)], in which ({r) represents the opening angle of 
the ejecta at radius r. So, equation(l) can be written as: 

— = 2Tm [1 - cos C(r)]r , (2) 

where E and Fq are the initial kinetic Energy and initial Lorentz factor of the ejecta {E = 
FqMoC^), and n denotes the number density of the cloud (p = nipn). Paczynski & Rhoads 
(1993) derived the relation between / and F from conservation of energy and momentum. 
Here, we use their relation in a form suitable for our computations: 



4f ^ _vnEi 

The ejecta's opening angle ({r) increases with increasing r as a result of lateral spreading 
of the cloud of ejecta+swcpt-up matter in the comoving frame at the sound speed c^, which 
has been derived by Rhoads (1998) to be as below: 

dCir) = (4) 
r 

where tco denotes the time from the event, measured in the ejecta comoving frame. Substi- 
tuting dtco — dt/r, and dt — dr/ (5c, we have : 

dC{r) ^ Cs/c 
dr pVr' ^ ' 
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Now, eliminating / between equations(2) and (3) yields: 



dr mpToc^. I v^r2 - 1 

— = — ZTrn— — 1 — cosC(r) — r . (6) 

Let's rewrite equations(5) and (6) in a non-dimensional form as below: 



dC, Cs/c 

dr, - r,r(,,)(l-J7)-i/2 

where, the non-dimensional parameter r\ is defined as: 



(7) 



in which: 

,(„/^) , = 4.75 X 10" (ll) cm . (9) 

These coupled first order differential equations can be solved numerically by introducing the 
initial conditions: 

where 770 = TQ/l{n/ E). 

Noting that (5 = dr/cdf = (1 - l/r2)V2^ we have: 



In equation(ll) we used equation(8) and a non-dimensional time parameter r defined as: 

ct , ^ 

' ^ WW) 

The numerical results of cquation(ll) are used in appendix C, where we consider the effect 
of burster geometry on the observable time duration. 



2.2. Formulation of Time Duration Distribution 



We begin with introducing the probability density for a collimated burst to occur in 
a direction through the cloud with a thickness L and a number density n. We assume the 
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cloud thickness to have a gaussian distribution in various directions from the central engine. 
By assuming a gaussian distribution for the cloud density as well we have: 



dn dL dD, Att v^ttu. 



exp 



— (n — ny 
2^^ 



exp 



(13) 



The quantities L and denote the mean thickness of the cloud and its dispersion, respec- 
tively. Likewise, n and are defined in a similar manner. Denoting the angle between the 
ejecta symmetry axis and the line of sight by 9 (Fig. 2), and considering the independence of 
the probability density from azimuth angle in equation(13), we can write: 



2ti sin e 1 



dn dL d9 



An 



2txo 



exp 



— n — n 



2a2 



exp 



(14) 



The synchrotron emission is a fast process for our model (see appendix A), so T.,.ec (the time 
duration of a GRB as measured by an observer cosmologicaly near to the source and located 
on the line of sight), is attributed to the time that the shock front takes to cross the dense 
cloud. As explained in appendix C, the cloud thickness L can be expressed as a function of 
n, and T^ec (Eqn.[C16]): 

L^L{Trec,n,e) . 



So we write: 



or: 



dn d9 



dp 
dL 



d^p 



d^ 



dp 



dn d9 \ dT, 



d^p 



dTr, 



dL 



dL 



(15) 
(16) 

(17) 

dn d9 d log Tree dn d9 dL \d log Tree J 

notifying that by this substitution(Eqn.[15]), the bursts that do not manage to cross through 
the cloud with a thickness L (and stop in it) are practically omitted (see appendix C). Using 
equation(14) in equation(17), it is seen that: 

d^p 



sin 9 1 



dn d9 d log T^, 



471" o-^ o". Vc^logT, 



dL 



X exp 



Integrating over 9 and n yields: 
dp 



dlogTr, 



1 r2 
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exp 
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(18) 



(19) 
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The effect of red shift is not considered yet. Equation (19) only gives the probabihty density 
for an observed burst to have a specified logarithm of time duration, as measured by an 
observer near to it. We now investigate the relation between dp / dAogTyec and dp/d\ogT^ 
where stands for the time duration of a GRB measured at Earth. To obtain the later, 
the former must be integrated over red shift z, using a weight function ^^^^(z), so that 
Fgrb{^) represents the probability for occurring a GRB in a red shift between z and 
z + dz. To show this, let's consider an observer located on the hne from us to an occurred 
GRB which is (cosmologicaly) near to it. The probability density for the GRB to occur in 
a red shift z (with respect to us) , and to have a specific log Tree (measured by the observer 
near to the GRB), is clearly as below: 



d'^p 



dz d log Tr, 



dp 



dlogTr, 



(20) 



To obtain d'^p/dz dlogT^, which is the probability density for observing a GRB occurred at 
a red shift z and observed to have a specific T®, we write: 



d'^p 



d 



dz d log dz \d log 



dp 



d_ 
dz 



dp 



dlogTr, 



dlogTrec 

dlogTe 



(21) 



noting that = (1 + z)Tree, the second term in the the bracket equals one. Now, using 
equation (20) we have : 



d'^p 



dz d log T^ 



dp 



dlogTr, 



(22) 



Trec=T(s/{l+z). 



After integrating equation(22) over z, the final form of the observable probability density 
will be as below: 



dp 



rflogTg: 



GRB 



dp 



z=0 



dlogTr, 



dz . 



(23) 



The exphcit form of -^^^^(-z) is needed. This form is obtained in appendix D (Eqn.[D4] and 
Eqn.[D5]). The second term in the integrand is given by equation(19), in which the implicit 
form of L{Trec,n,6) appears. Appendix C is devoted to the procedure of obtaining this 
function. For evaluation of dp/dlogT^ (Eqn.[23]), we need the values of eight parameters. 
Four of them are the cloud parameters L, a^, n, and (T„, and the fifth one is the index q 
corresponding to the GRB occurrence rate (see Eqn.[D5]). The later three ones are the ejecta 
parameters, E, Fq, and Co, which are the initial kinetic energy, the initial Lorentz factor, and 
the initial opening angle of the ejecta, respectively. 
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NUMERICAL COMPUTATIONS AND RESULTS 



The Mathmatica 4 software is used in the numerical computations. In the procedure 
we begin with solving the coupled differential equations(7) and (If) which govern the ejecta 
evolution, and in which Fq and Co are the only free parameters. Furthermore, we take 
r)o = 0. For fixed values of these parameters, the functions F(?7), r(r/) and Cradiv) (see 
Eqn.[C3]) can be uniquely obtained . Obviously F decreases with increasing rj (Fig.3), and 
reaches to F = 1 when r) approaches a certain value. This in fact takes an infinite time 
and results in an infinite non-dimensional time duration Tree (see Eqn.[C4] and Eqn.[C5]). 
So, to circumvent the difficulty, we consider an effective lower limit Tmin for the Lorentz 
factor of shocked matter, which yields an effective upper limit for non-dimensional time 
r. We adopted F^m = 2 as a lower cut-off. Then, following the procedure explained in 
appendix C, for a specific cloud thickness L, and correspondingly a specific i]^ (Eqn.[C13]), 
the non-dimensional time duration Trec{&,VL) ^ GRB can be calculated for every 9 and 
T)^ (Fig.4). Let's denote the radius corresponding to F^j„ by rjm so that F^j„ = r{r)jn), and 
recall that in the model, for r)m < Vlj emitted photons would be completely scattered by 
the electrons of not-swept part of the cloud (that they have to pass through, before entering 
free space; see appendix C), and so, the whole phenomenon may be called a "failed GRB" . 
But, if r]rn > Vl' shocked matter succeeds to go out of the cloud and, as explained in 
§ 4.3, due to a suppression process that intensively decreases the cross-section of Compton 
scattering, (the major part of)the emitted photons finally succeed to get released from the 
shocked medium and enter free space, provided that C(^l) > [\/V3 F(?7^)]^^ (see appendix 
B). It is really for this reason that Tree happens to be a function of 77^, and practically 
independent of 77^. Then, solving Tred^jVL) Vl numerically, we can obtain the function 
r]^{Trec,d) which is the equivalent non-dimensional form of expression (15). We rewrite 
equation(19) in a non-dimensional form suitable for numerical computations: 



dp 



d log Tree '^Q<iJe=0 



exp 



i/=0 



2^2 



exp 



X 



l{nv/E) ( df] 



dlogTr, 



2<j2 



sin 9 d9 du, 



(24) 



in which, equation(C13) and the definitions below: 
u = n/n , 
g = an/n , 
and 

<; = gl/L 

are used. So, after choosing the quantities L, cr^, n, (7„, and of course we can evaluate 
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the integral appearing in equation(24). Then, after choosing a value for q, the integral of 
equation(23) can be performed to obtain the observable quantity dp/dlogT^. 
As seen in equation(24), the initial kinetic energy E and the cloud mean density n appear 
only in the form E/n, and therefore they can not be found separately in a fitting process, 
and all that can be obtained is only their ratio. But, the observed fiuence of GRBs reveals 
that the released energy in GRBs is of order of Eiso ~ lO^^ergs for a isotropic burst, which 
reduces to {il/A7r)Eiso if the bursts were confined to a cone with a solid angle Q. We used this 
amount of isotropic energy to relate E to Co {E = EisoQI'^), and reduce the free parameters 
of the model to seven as Fq, Co, a^,n, an and q, which hereafter are called U parameters. 
These parameters must be chosen so that the results make the best fitting to the observed 
distribution of GRBs. To achieve the task, one must search in the seven dimensional 11 
space, and find the point in which the statistical quantity takes the smallest value Xmin- 
We used the BATSE 4th catalogue (Paciesas ct al. 1999) of 1234 GRBs and adopted the 
bins AlogT0 = 0.2 in a range —1.9 < logT® < 2.9, and used the gradient search technique 
to move toward the best point in which gets minimized. The obtained fitted values are: 

To = 0.97 X 10^ Co = 0.01 

L = 1.7 X lO^^cm n = 2.9 X lO^^cm"^ 

= 0.21L an = 0.71ri 

q = -0.70 (25) 

with a corresponding value Xmin — ^-^ (pcr degree of freedom). As can be seen in Fig.5, the 
deviation from a more perfect coincidence to data seems for the structure in the observed dis- 
tribution located around logT® ~ 0.7. The structure has been noted before (Yu et al. 1998). 
We put aside the the data of the noted structure, which are ones with durations between 
0.3 < log(T'0) < 0.9, and again repeated the numerical searching process in the parametric 
space. The new obtained values of the fitted parameters are: 



To = 0.96 X 10^ Co = 0.01 

L = 1.8 X lO^^cm n = 2.9 X lO^^cm"^ 

a^ = 0.211 an = 0.7m 

q = -0.70 (26) 

which slightly differ from what previously obtained (Eqn.[25]). But this time, Xmin reduces 
to 1.1 (Fig. 6). So the mentioned structure may be interpreted as a result of an independent 
phenomenon or effect which was not considered in our model. 
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4. DISCUSSION 

4.1. The GRB Source 

Though the original shapes of equations (13) and (14) are naturally normalized, the 
resulting final equations (19) and (23) are not, because: 

1) the bursts that their symmetry axes make an angle 6 > CradiVL) ^e detected (see 
appendix C), 

2) thc bursts for which 77^ < f]^ were not considered in the numerical computations (because 

they do not manage to cross out the cloud). 

The total probability of observing the occurred bursts, j'_^{dp / dlogTco) dlogT®, is obtained 
to be 1.47 X 10~^, using our best fit parameters(Eqn.[26]). This small probabifity must be 
interpreted to be due to the above two reasons. The first reason describes the suppression of 
observed GRBs by the term (Q/47r) ~ Q/2 ~ 10~^, while the second is responsible for the 
remaining factor of 10~^. So, only about one percent of the bursts manage to produce a real 
GRB, and only about 10~^ of these GRBs occur in our line of sight. Of course, the lateral 
spreading of the ejecta+swept mass may modify these two factors, increasing the first and 
decreasing the second. 

With the values for the fitted parameters in equation(25) or in equation(26), the mean mass 
of the clouds M = ^TrnnipL is about 1.2 x 10^'' gr ~ GMq, which is of the order of the 
envelop mass in massive stars. Such amorphous clouds seem strange in stars, but in close 
neutron star-supergiant binaries where the neutron star orbits around the core and accrete 
the envelope, the spherical symmetry is likely removed, as pointed out by Podsiadlowski et 
al. (1995). Terman, Taam, & Hernquist (1995) show that the system would emerge to form 
a red supergiant with a massive Thorne-Zytkow Object (TZO) (Thorne & Zytkow 1977). 
Podsiadlowski et al. (1995) also estimated a TZO birth-rate of > 10^'^yr^^ in the galaxy. 
Considering equation(D4) and the obtained total probability of observing a "real" GRB, 
which is 1.47 x 10^^ in our model (see above), it is seen that the GRB observation rate (2 
events per day) implies a total (real+failed) GRB rate of the order of ~ 10^^ Mpc^^yr^^. 
This is not too far from what Podsiadlowski et al. (1995) theoretically estimated for TZO 
birth-rate (> IQ'^galaxy'^yr'^). Qin et al. (1998) introduced AICNS (Accretion-Induced 
Collapse of Neutron Stars) scenario as GRB engines. Katz (1994) introduces a dense cloud 
model to explain the observed Gev gamma-rays in a number of GRBs (Dingus et al. 1993; 
Jones et al. 1996, e.g.), and suggests the amorphous envelops of TZOs for these clouds. 
The initial opening angle of the ejecta in our model (Co — 10~^ra(i ~ 0.6°) is obtained 
during the fitting process (§ 3). Such a small opening angle for an ejecta or a jet might 
be explained by attributing it to the collimating process of an ultra-relativistic ejecta with 
Fo > 1/Co in a sufficiently high magnetic field (Begelman,Blandford,&Rees 1984). The initial 
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Lorentz factor of the ejecta, as obtained in our model (Fq ~ 10^, see Eqn.[25]) provides the 
necessary condition of Fq > 1/Co- Aside from these theoretical justifications for the idea of 
a highly collimated ejecta at the source, there are some pieces of evidences supporting this 
idea (Lamb,Donaghy,& Graziani 2004; Waxman 2003; Granot 2003). 

4.2. The Bimodality 

At the mean time the observed bimodality must be interpreted as a result of the second 

reason expressed in § 4.1. Though in our model one may expect only one heap in the logT^ 
distribution associated with the directions in the clouds having both the most probable L 
and n (which are equal or near to L and n), but as the numerical computation shows, such 
directions in the clouds arc too thick and too dense to be crossed out by the ejecta, and 
therefore a real GRB would not be produced. So, we are left with four other regions with 
high probabilities: 

1) n ~ ri and L < L (short GRBs) 

2) n ~ n and L > L (no GRBs) 

3) n >n and L L (no GRBs) 
A) n <n and L ~ L (long GRBs) 

Now, we show that the directions trough the clouds associated with region(l) produce the 
short duration GRBs, and ones associated to regions (2) and (3) produce no GRB, while the 
others in the forth region produce the long duration ones. 

As to the equation(9) the quantity l{n/E) (which is of the order of the sedov length) is 
~ 7 X 10^-^ cm <^ L for n ~ n (see Eqn.[25]). As seen in Fig. 7, the time duration of GRBs 
associated to such directions is of the order of the time duration of short GRBs (Case (1) 
above). Fig. 8 shows that in our model the calculated time duration of GRBs for n ~ lO-^ri 
is of the order of the duration of long GRBs. Using equation(9), we see that in this case ( 
n ~ 10~^n ), l{n/E) is about 5 x lO^^cm ~ L. So in our model the long duration GRBs 
are due to the passing of ejecta through the directions where n <C n and L ~ L (case (4) 
above). Furthermore, Since l{n/E) <^ L when n ~ n, we see that in cases (2) and (3) 
above, the ejecta would stop in the dense cloud and the produced photons can not scape 
from the optically thick cloud. In Fig.9, Trec{L]n,9) is plotted for a number of densities, 
when 77^ = L/l{n/E) has the highest permitted value 77^. 

These general features of our calculations result from the general features of our model, and 
therefore we speculate that any distributions for the clouds thickness and density which are 
picked around a mean value could explain the general features of the duration distribution. 
Our chose of gaussian distributions for thickness and density was only for the few parameters 
needed to describe them. 
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4.3. The Opacity 

In "The Dense Cloud Model" at this stage we have simply omitted the bursts which 
ejectas can not go out of the cloud and stop in it (because the emitted photons would be 
scattered by the dense cloud) , and we have claimed that the produced photons by all bursts 
that succeed to cross out the dense cloud can finally enter the free space. 

At the first glance the model might appear to have a serious problem in opacity, namely, 
in this model the best fit mean density and mean thickness of the clouds are found to be 
of the orders of n ~ lO^^cm"^ and L ~ lO^^cm. So, as to the relation Top = (Tj, n L, 
= 6.65 X one would expect an optical depth of the order of 10^. But there 

are two factors that remedy the situation: 

(i) As the numerical calculation shows (see § 4.2), the most probable directions character- 
ized by n ~ ri and L ~ L are dynamically too thick to be crossed out by the ejecta and the 
radiation produced in this case would be completely scattered by the dense cloud. On the 
other hand, as discussed in § 4.2, the long duration GRBs are due to the crossing of ejecta 
through directions where n ~ IQ-^n and L ~ L. since the density is reduced by the factor 
10^^, the optical depth drops to the order of 1. As explained in § 4.2, the short duration 
GRBs are due to the directions in the cloud where n ~ n and L ~ 10~^L. In this case the 
optical depth of the cloud reduces to ~ 10"^, which is yet too high. But, 

(ii) most of photons emitted off the shock front have the chance to be overtaken by the mov- 
ing shock (appendix B). Moreover, it has been shown (Shekh-Momeni & Samimi 2004) that 
in high temperatures /cT ~ 10^ mgC^ the cross section of Compton scattering for instance 
for Mev photons effectively drops to ~ 10~^ a^, so such a high temperature plasma is much 
more transparent than what may seem at first. The same is true for the case of a power 
law distribution of electrons (see Eqn.[A4]). Most of the electrons in such a distribution 
have energies of the order of 7e.min^eC^ ~ 2 x lO^nicC^j in which we have used equation(A5), 
taking p = 2.5, .^e = 1/3 and F = 10^. For this case the Compton cross section suppresses by 
the factor of ~ 10'^ — 10~^ for lOOet' — IKev photons (the energy of the emitted photons in 
the shocked medium comoving frame is less than their observed energy by the factor of F), 
and therefore, the optical depth drops to the order of 1. So the photons overtaken by the 
shocked medium may remain in it without being scattered until the shocked medium crosses 
up the dense cloud. Briefly, due to this second factor, the optical depth in short duration 
case diminishes to Top ~ 1, and the optical depth in long duration case diminishes to Top <^ 1. 
The variability seen in the light curves of long duration CRBs maybe attributed to the het- 
erogeny of cloud's density in the ejecta's trajectory. This maybe the case in long duration 
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GRBs for which Top <^ 1, but in short duration GRBs for which Top ~ 1, we expect that such 
a heterogeny would not be appeared. 

The lateral expansion of the shocked medium can reduce its density and therefore its optical 
depth. This effect of course may cool the shocked medium and cause an increase in effective 
Compton cross section in it. The consideration of dynamical feedback makes our modeling 
much more complicated and it has been neglected here. In "The Dense Cloud Model" at 
this stage we have simply omitted the bursts which ejectas can not go out of the cloud and 
stop in it, but have claimed that the produced photons by all bursts that succeed to cross 
out the dense cloud can finally enter free space. 

A complete study needs to include flux computations. An exact comparison with ob- 
served time duration data would be possible only when the theoretical time duration ap- 
pearing in this work were exactly evaluated in the same manner as the observable dura- 
tion Tgo is defined (as the time interval over which 5 percent to 95 percent of the burst 
counts accumulate). Moreover, the BATSE's triggering mechanism made it less sensitive to 
short GRBs than to long ones and therefore short GRBs were detected to smaller distances 
(Mao,Narayan,& Piran 1994; Cohen & Piran 1995; Katz & Canel 1996), so a smaller num- 
ber of them have been observed. Lee & Petrosian (1996) studied this effect and corrected 
the number of short GRBs. In a more exact work this correction must be considered too. 

We acknowledge the anonymous referee for valuable comments. F. S. acknowledges Mr. 
Mehdi Haghighi for his guides on computational methods. This research has been partly 
supported by Grant No. NRCI 1853 of National Research Council of Islamic Republic of 
Iran. 
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A. Scyncrotron Cooling Time in Dense Cloud Model 

In this part we use the review paper of Piran (1999) to estimate the synchrotron coohng 
time in our modeL The synchrotron coohng time in the comoving frame is: 

tsyn,co — . J J 1 (Al) 

where ut is the Thompson cross section and 7e is the Lorentz factor of the emitting electron, 
while Ug stands for the energy density of magnetic field which in GRB literature is assumed 
to be proportional to the comoving internal energy density of the shocked matter: 

= , (A2) 

so that represents the share of magnetic field in m, which is given by: 

u = Ar^nrUpC^ , (A3) 

in which, n is the number density of surrounding medium (ISM, or a dense cloud as assumed 

in our model), and F denotes the Lorentz factor of the shocked matter. 

The electrons in shocked media are assumed to develop a power law distribution of Lorentz 

factors: 

A^(7e) CX for : 7e > 7e,mm • (A4) 

The convergence of total energy of the electrons requires the power index p to be greater 
than 2, while the assumed lower limit 7e,min is to prevent the divergence of electron number 
density and is obtained to be: 

_rnp p-2 
me p - 1 

where = Ue/u represents the share of the electrons in the internal energy of shocked 
matter. 

Furthermore, the time interval between emission of two photons from the same point in the 
comoving frame, 5tco, and the time interval 5t0, which represents the time interval of their 
successive arrival to a cosmologically distant observer (at the earth), are related as below: 

<^^e = (l + ^)%- (A6) 

Equations(Al-A6) are adopted from Piran (1999). Now, by substituting equations(A2),(A3), 
& (A5) in equation(Al), and then using equation(A6), the synchrotron cooling time mea- 
sured by a terrestrial observer turns out to be: 

< |(i + .)C'G' (^)' (^) (..«.)-'r-.ec 

~io-(i + .)(|^)"(f|)"(^)"'r-W, (A7) 
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which is clearly much less than all observed time durations of GRBs. So, the synchrotron 
emission in a dense cloud model must be considered as an instantaneous process, and there- 
fore the shock front must be regarded as the emitting surface. This allows us to attribute 
time duration of a GRB merely to the time that the shock front takes to cross the dense cloud. 



B. The relation between Lorentz factors of the shocked matter and the 

emitting surface 

Here, we want to find the relation between F and T' which are respectively the Lorentz 
factors of the shocked matter and of the shock front (which is the emitting surface in our 
model; sec appendix A). In Fig. 10, j3 and j3' correspond to shocked matter and shock front 
speeds, both measured in the (central) source frame. In the shocked matter frame (Fig. 11) 
the dense cloud which is seen to have a density F n, moves toward the shocked medium with 
a speed /3, while being compressed to a density equal to 4Fn. Consequently, as illustrated in 
Fig.ll, the shocked medium expands towards right with a speed P'^^ which is in fact the speed 
of the shock front in the shocked matter frame (compare Fig. 10 and Fig.ll). Considering 
the conservation of nucleon number, it is seen that: 

Peo = ^- (Bl) 

Considering the relativistic summation of velocities, we have: 

Noting the relations F = (l-/32)-V2 and F' = (l-/3'^)-V2^ the substitution of equation(Bl) 
in equation(B2) finally yields: 

The distance from the origin to the shock front, denoted by r', can be obtained by integrating 
(3' over t: 

r' = I (5' cdt + ro . (B4) 
Jo 

The results are shown in Fig. 12. As seen, the difference between r and r' never exceeds one 
percent. This result justifies the use of symbol r (instead of r' )for the location of the shock 
front through out our formulation. 

Now, let's consider a photon that is radiated from the shock front (the emitting surface) 
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moving with the Lorentz factor F' (Eqn.[B3]), in a direction which makes an angle e with 
the velocity vector of the shocked matter (Fig. 10). if : 

e > T'-\ (B5) 

the photon would be overtaken by the shock front. The majority of the emitted photons 
fulfill this condition when F' ~ 10^ and Co ~ 10"^ (see Eqn.[25]). 



C. Geometrical Considerations 

Here the effect of burster geometry on the observed time duration is investigated. At 
first, as shown in Fig. 13, we consider a radiating segment on the shock front. The symmetry 
axis is denoted by z'. Noting the definition of e in the figure, it is seen that in addition to a 
radial component dr/dt, the velocity vector of the segment must have a lateral component 
viat which is equal to: 

as measured in the source frame. In equation(Cl) the non-dimensional radius rj, as defined 
in equation(8), is used instead of r, . Equation(Cl) is obtained simply by using equation(5) 
and assuming that the lateral speed of the segment in a frame moving (only radially and) 
instantaneously along with the segment, is the fraction e/({r]) of sound speed Cg (which is the 
lateral speed of the emitting surface at its edges), and noting that the lateral speed in the 
source frame is less than its corresponding value in the (radially instantaneous) comoving 
frame by the factor 1/F. It is well known that the radiation emitted by the segment is almost 
confined to an angle 1/F. In Fig. 13 the axis of the radiation cone emitted by the segment 
is denoted by z", which is parallel to the velocity vector of the segment and, as seen in the 
figure, makes an angle e + 5 with z' axis, where 5 — tan'^leCg/lrpcQ]. So the radiation 
angle Cradiv^ ^) depicted in the figure) can be written as below: 

C„.to, e) = e + tan- + Um . (C2) 

and the total radiation angle, defined as the radiation angle at the edge of the emitting 
surface, can be written as: 

This angle represents the cone of space illuminated by the emitting surface. Using the results 
of equation(7), the total radiation angle can be evaluated for every "77". Having defined the 



-16- 



total radiation angle, we explore the effect of burster geometry on its observed time duration 
Tq. As shown in Fig. 2, the problem is studied in a spherical coordinate system in which the 
central engine is taken as the origin, and the line-of-sight as the polar axis z. Consider a 
photon emitted off a point on the emitting surface (the shock front) with radial coordinate 
r and polar coordinate 0; and at an instance t, measured in the source frame (the point 
is not shown in the figure). The relation between t and the photon arrival time tree to a 
(cosmologically) near observer is obtained by Granot, Piran, & Sari (1999). Making suitable 
for our model, it is adjusted to the form below: 

t,„ = t_!li5^£|^. (C4) 

In this equation the instance t = is defined as the time that the ejecta collides with the 
dense cloud at r = ro; while tree = is the time that the (cosmologically) near observer 
receives the photon emitted ai t — from the point with coordinates r — and © = 0. 
Defining: 

Ctrer ^^^^ 



l{n/E) ' 

and noting equation(8) and equation(12), we rewrite equation(C4) in the form below: 

Tree — T — T) COS Q + TJo . (C6) 

Multiplying cquation(C4) by the cosmological time dilation term (1 + ^), results in the arrival 
time as maybe observed at the earth: 

V ^ ( rcosO — ro\ ,^„, 

= (1 + Z%ee = (1 + -2) U ^ ^ j , (C7) 

or equivalently: 

T© = (1 + Z)Trec = (1 + 2;) (t - 7? COS © + 770) , (C8) 

in which we used the non-dimensional time duration r© defined as below: 

"^^^ (C9) 



l{nlE) 



Now, as shown in Fig. 2, we consider a situation where the ejecta's symmetry axes makes an 
angle Q with the line of sight. The necessary condition that at least some photons of the 
emitting surface are detected by the near observer is: 

Uiij]) > e . (CIO) 



Here we denote the inverse of function Cradiv) by rjradiCrad) , which gives the radius corre- 
sponding to (rod- As seen in Fig.2, for ^'s larger than Cradivo): the first photons reaching the 
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detectors are those emitted at 77 = r]rad{d)- So, we use T{r)) (which is obtainable by solving 
equation(ll)) to write: 



if 9 < CradiVo) 
r{r]rad{0)) if 6* > CradiVo) 



(Cll) 



where ti{9) represents the starting time (in the source frame) that the emitted photons can 
reach the (cosmologicaly near) observer. Now, using the numerical results of equation(ll) 
we can obtain the function r]{T). Then, considering equation(C6), the non-dimensional time 
Trec,i{9) corresponding to ti{9) would be as below: 





r{l]rad{9)) - r}rad{9) COs[9 - C(^rad(6'))] + Tjo 



0<Co 

Co<9 < CradiVo) 
9 > CradiVo) 



(C12) 



Now, we are to find the time Trec,2{9) after which no photons can be detected by the near 
observer. In Fig. 2, the photons emitted from the edge point A can reach us at all times 
greater than ti{9). Let's remind that in our model the emission process terminates at the 
time that the shocked matter goes out of the cloud (appendix A). As is seen in equation(C6), 
the closer to the point S is a point on the emitting surface, the later its emitted photons 
would reach the near observer, of course, provided that the observer line-of-sight remains in 
the radiation cone of the emitting point. 
Defining: 



L 



Vl = 



l{n/E) ' 

and recalling equation(C2), we can solve the equation: 



(C13) 



(C14) 



to find the function e = e(r7^, 9)^ which specifies the furthest point (on the shock front at the 
radius 77^) which its radiation reaches us, of course, if it remains smaller than C(^l)- 
Now, using equation(C6), we can find the instance that the last photons reach the near 
observer: 

rrec,2i0,VL,Vo) ^ r{r]^) - T]^ cos{9 + min[e{r]^,9),CiVL)]) +Vo ■ (C15) 

Finally, the non-dimensional time duration of a GRB, Trec{9, r/^), will be equal to Trec,2{9, rj^) — 
Trec,i{(^) and, as to equation(C12) and equation(C15), besides being a function of parameters 
of the ejecta and the cloud, it is also a function of the inchnation angle 9. So, the time 
duration of a GRB as measured by an observer cosmologivally near to it would be a function 
of L, n, and 9: 



Tree — T^ee^L, 71, 9^ , 



(C16) 
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where: 

Trec^TreMv,)^-^^^, (C17) 

The expression (CI 6) needs some explanations. In our model the ejecta parameters Co : E 
and To and the radius tq are assumed to be the same in all GRBs. So, these parameters do 
not appear in equation(C16) explicitly , though it is implicitly a function of them too. At the 
mean time, the cloud's density n and its thickness L are assumed to be different in different 
directions, and therefore they do appear explicitly in the expression (C16). Furthermore, if 
a cloud thickness L is much more than its associated Sedov length Isedov = {E/niiip c^)^/^, 
the ejecta may not cross it up, and finally stops in it. In such a situation the time duration 
of a GRB would not be a function of L, while it still remains dependent on n and 9. This 
is why in expression (C16), the quantity L is distinguished from n and 6' by a semicolon. 
The rearrangement of the expression (CI 6) in the form L = L{Trec, )0, n would be meaningful 
only when we are dealt with the situation where the ejecta succeed to cross up the cloud 
(see Eqn.[15]). 



D. The explicit form of F^j^g{z) 
Here the relation between F^^^gi^z) appearing in equation(23) and the GRB occurring 



GRB 

rate fajigiz) (in units of Mpc~^ yr~^) is derived, so that by adopting a cosmological model 
for the occurring rate, the integral in equation(23) can be evaluated. 

The number of GRBs that their effects could reach us in a time interval Sto (which is very 
much less than the present comoving time to), and from a spherical volume clement SV^ 
(disregarding various effects, such as the geometrical ones discussed in appendix C, or those 
related to detectors threshold which affect the number of detected GRBs) is as below: 

i- + z 



St, 



f,,,{z).inR'{z)r{z)Sr{z)-^, (Dl) 

i. + z 

where r{z) denotes the non-dimensional radial parameter of the source that its effects reach 
us with a red shift z, and R{z), the scale factor at this red shift. In Einstein-dc Sitter model 
we have: 

i?or(z) = -f{l-(l + z)-V2}, (D2) 
where Hq is the Hubble constant. Furthermore, in FRW metrics: 

R{z) = Ro{l + z)-^ , (D3) 
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where Rq = R{z — 0). Now, using equations(D2) and (D3) in equation(Dl) we obtain: 

3 

Fansiz) oc = (1 + z)-'^' {1 - (1 + zyy^y ■ m 

What is remained is the exphcit form of fangiz). The high variabihty seen in GRB hght 
curves has convinced the investigators to relate GRBs to stellar objects , and consequently 
their rate fa^giz) to the star formation rate fgp{z). The simplest model is of course a 
proportional model f^j^giz) oc fgp{z). the proportional model may be correct if GRBs are 
attributed to the evolution of massive stars whose lifetime is negligible in comparison with 
the cosmological time scale, but in NS-NS mergers model the proportionality may not be 
valid (because of the delay time from the star formation to NS-NS mergers). Wijers et 
al. (1998) claimed that there is a good consistency between the proportional model and 
the observed GRB brightness distribution, while Petrosian & Lloyd (1998) concluded that 
none of the NS-NS and the proportional model are in agreement with the observed fgp[z). 
Totani (1999) ascribed this discrepancy to the uncertainties in SFR observations. Anyway 
we simply assume the GRB rate to be as below: 

/chb(^)=/c«b(0)(1 + ^)'^'' , (D5) 

and treat g as a free parameter that its best value should be obtained during the fitting 
procedure. In Fig. 14, Fq^j^{z) is plotted for a number of g's. Clearly the case g = 
corresponds to a universe where the changes in the rates of astrophysical phenomena are 
only due to its expansion (non-evolutionary universe). It can be seen that -Fg^b(^) takes its 
maximum at z ~ 1, which is not very sensitive to the magnitude of q. 
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Fig. 1. — Distribution of log of GRBs' time Duration (normalized), obtained from data of 
BATSE 4th catalog. 
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Fig. 2. — Geometry of Radiation. In the figure, r]rad{(^) denotes the radius where the radiation 
angle (rad becomes equal to 9 (see appendix C). As shown in the figure (for 6''s larger than 
CradiVo)) the first photons reaching us are those emitted from the edge of the shock front at 
this radius. If e{r)j^,9) < C(^i,)) then the last photons reaching us are those emitted from 
e = e{rij^,9) on the shock front (point S in the figure). The point S would be the furthest 
visible point on the shock front, since the line of sight fells out of the radiation cone of the 
points farther than it (points between S and B). If e(?7^, 9) > C(^i)) then the point B would 
be the last visible point (see Eqn.[C15]). 
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Fig. 3. — The Lorentz factor of shocked medium r(77) versus the non-dimensional radius rj, 
with To = 1000 and Co = 0.01 (see Eqn.[7]). 




Fig. 4. — The cvahiatcd non-dimensional time duration Trec{(^, Vl) GRBs, versus 6 and 7]^, 
with Fq = 1000 and Co = 0.01, (see appendix C). 
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Fig. 5. — Results of best fitting. The sofid curve shows the observed time duration distri- 
bution and the dashed one is calculated using the best of values for the model parameters. 
xLn = l-W(see§3). 




Fig. 6. — Results of best fitting. Same as Fig. 5 with the data corresponding to the structure 
with 0.3 < logT© < 0.9 were cut. = l.lpdf (see § 3). 
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Fig. 7. — The time duration of GRBs, evaluated in our model with n — n. It is plotted 
versus the inclination angle 6 for a number of possible thicknesses 77^. As seen, the order of 
time durations are of the order of ones in short duration GRBs (see § 4.2). 




Fig. 8. — The time duration of GRBs, evaluated in our model with n = lO^^n. It is plotted 
versus the inclination angle 9 for a number of possible thicknesses 77^ . As seen, the order of 
time durations are of the order of ones in long duration GRBs (see § 4.2). 
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Fig. 9. — \ogTrec{L; n, 6) is plotted versus the inclination angle 6 for a number of densities, 
with = L/l{n/E) equal to the highest permitted value rj^- The number near each curve 
is log(n/n) (n = 2.9 x 10^'' cm~^) (see § 3 and § 4.2 ). 



swept-up mass cloud 




(emiting surface) 



Fig. 10. — The shocked matter and the shock front (emitting surface) speeds, as viewed in 
the source frame. In the figure, if the angle e were larger than l/F' the photon would be 
overtaken by the shock front (from which it were emitted) , so that the photon would remain 
in the shocked medium (see appendix B). 
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swept-up mass cloud 



CO 



shock front 



Fig. 11. — The shock front (emitting surface) speed measured in the shocked frame. In 
this frame the cloud has a density of Tn and moves toward the shocked medium with a speed 
P, and becomes compressed to a density of AFn. So clearly we have /S^^ — P/^ (see appendix 
B). 




Fig. 12. — (r' — r)/r versus 77 = r/l{n/E), with n = 2.9 x lO^^cm ^ (see appendix B). 
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Fig. 13. — Geometry of Radiation. The lateral speed of the emitting segment, viat, causes 
the radiation cone axis z" to make an angle e + 5 with z' axis, where 5 — tan~^[eCs/ (r/3cC)] 
(see appendix C). 
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Fig. 14. — The quantity -Fghs(-^) (the probability density of observing a GRB at a redshift z) 
is plotted for a number of g's (the value of q is written near to the peak of its corresponding 
curve). As seen, all curves have a maximum at 2; ~ 1. This explains why the most of GRBs 
have redshifts z ^ 1 (see Eqn.[D5]). 



